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We compute the transition temperature Tc and the Ginzburg temperature Tq above Tc near a 
quantum critical point at the boundary of an ordered phase with a broken discrete symmetry in 
a two-dimensional metallic electron system. Our calculation is based on a renormalization group 
analysis of the Hertz action with a scalar order parameter. We provide analytic expressions for 
Tc and Tq as a function of the non-thermal control parameter for the quantum phase transition, 
including logarithmic corrections. The Ginzburg regime between Tc and Tq occupies a sizable part 
of the phase diagram. 
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I. INTRODUCTION 

Instabilities of the normal metallic state lead to a 
rich variety of quantum phase transitions in interacting 
electron systems. Near a quantum critical point elec- 
tronic excitations are strongly scattered by order param- 
eter fluctuations such that Fermi liquid theory breaks 
downiii^ The fluctuation effects and the ensuing non- 
Fermi liquid behavior is particularly pronounced in two- 
dimensional systems. It is therefore not surprising that 
quantum critical fluctuations are frequently invoked as 
a mechanism for the enigmatic strange metal behavior 
observed in cuprate superconductors and other layered 
correlated electron compounds. 

Theoretical works have mostly focused on the quan- 
tum critical point and its extension into the quantum 
critical regime at finite temperature, where quantum fluc- 
tuations are particularly important. Less attention has 
been paid to the Ginzburg region near the critical tem- 
perature, which is characterized by strongly interacting 
classical order parameter fluctuations)^ This is somewhat 
unwarranted since classical critical fluctuations also affect 
electronic excitations very strongly. In two dimensions 
they lead to a contribution of the order to the quasi- 
particle decay rate, where ^ is the diverging correlation 
length 4"— 

In this paper we compute the size of the Ginzburg 
region above the critical temperature near a quantum 
critical point in two-dimensional metals. More specifi- 
cally, we consider continuous quantum phase transitions 
associated with the spontaneous breaking of a discrete 
symmetry, described by an effective Hertz action^ for a 
scalar order parameter with dynamical exponents z = 2 
or z = 3. We calculate the transition temperature 
as a function of the non-thermal control parameter for 
the quantum phase transition, as well as the Ginzburg 
temperature Tq above T,.. The size of the Ginzburg re- 
gion Tq — Tc is determined to leading order in the dis- 
tance from the quantum critical point. The dependence 



of the Ginzburg temperature on the control parameter 
was derived already by Millis.- However, that study did 
not access the Ginzburg region between Tq and Tc- A 
comprehensive analysis of all finite temperature transi- 
tion and crossover lines in dimensions d > 2 was per- 
formed by Sachdevii^ For discrete symmetry breaking 
in two dimensions the critical temperature Tc and the 
Ginzburg temperature below Tc were recently compared 
within a renormalization group study which allowed to 
approach the finite temperature transition. However, 
in that work the flow equations were solved only numer- 
ically, while we now present analytic results. 

The paper is organized as follows. In Section II we de- 
rive the renormalization group equations for the effective 
Hertz action. These are solved analytically in an approx- 
imate form in Section HI. In Section IV we discuss the 
results for Tc, Tq and the size of the Ginzburg region, 
before concluding in Section V. 



II. HERTZ ACTION AND FLOW EQUATIONS 

Our analysis is based on the Hertz action^ 







(1) 



where 0(x, r) is a real scalar order parameter field and 
4'<i,u}n its momentum representation; ujn — 2'KnT with in- 
teger n denotes the bosonic Matsubara frequencies. For 
the dynamical exponent z we consider the cases z — 2, 
which describes density wave transitions, and z — rele- 
vant for a nematic transition or Ising-type ferromagnetic 
transitions. We do not address the issue under which cir- 
cumstances the Hertz action provides a faithful descrip- 
tion of quantum criticality in two-dimensional metalsJ^ 
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Before embarking on the renormalization group ap- 
proach, we would hke to emphasize that in two dimen- 
sions Tc{So) cannot be obtained from a first order ex- 
pansion in the quartic couphng uq, even if it is weak and 
irrelevant at the quantum critical point. To leading order 
in uq, the inverse susceptibility S is given by 



= Oq + ai y / 7-— :t 



|q|-2 



(2) 



where a is a positive constant. At finite temperature 
the Matsubara frequencies are discrete and the classical 
fluctuation contribution from a;„ = diverges logarith- 
mically in the limit (5o ^ in two dimensions. Trying to 
treat this divergence by a self-consistent equation, replac- 
ing So by S under the integral, one finds that the transi- 
tion temperature Tc is suppressed to zero at the critical 
point given hy 6 — 0, irrespective of Sq. This behavior 
is reminiscent of the Mermin- Wagner theorem, which ex- 
cludes spontaneous breaking of a continuous symmetry 
in two dimensions. However, the above first order calcu- 
lation is essentially independent of the symmetry of the 
order parameter, and is therefore misleading at least in 
the case of a discrete symmetry. 

We solve the problem by using flow equations which 
describe the renormalization of the inverse susceptibility 
(or "mass") S and the quartic coupling u due to fluctu- 
ations. The flow equations are derived from an approxi- 
mate ansatz for the exact effective action r'^[0], that is, 
the generating functional for vertex functions in the pres- 
ence of an infrared cutoff A.^^^ The cutoff is implemented 
by adding a regulator term of the form i J (pR^cf) to the 
bare action iS[0]. The exact flow of T^[(t)] is given by the 
Wetterich equationi^ 



dAT^[<j>] = -tr 



dAR' 



1 



2 r(2)AUl+i?A 



(3) 



where r'^^^'^[0] is the matrix of second derivatives of r'^[0] 
with respect to (p and the trace sums over momenta and 
frequencies. We approximate r'^[0] by an ansatz of the 
form Eq. ([T]) with a renormalized mass term 5^ and a 
renormalized coupling u^. Inserting this ansatz in the 
exact flow equation for r'^[0] and comparing coefficients, 
one obtainsi^ 



aAi?^(q) 



Saw'' = 3(u^)'r 



(27r)' 



5^ 



dKR^'ici) 



i?A(q) 



(4) 



i?A(q) 



(5) 

The initial conditions for the fiow are 6^° — 5q and — 
Uq, where Aq is a (fixed) ultraviolet cutoff. As a regulator 



we choose the Litimi^ function -R'^(q) — (A^ 
q2), with derivative 9Ai?'^(q) = 2Ae(A2 



-q2)e(A2- 
q^), which 



restricts the momentum integrals in the flow equations 
to |q| < A and replaces the q^-term in the denominators 
by A^. From now on we fix the dimensionality to d=2. 

The Matsubara sums in the above flow equations can 
be expressed in terms of polygamma functions ^'„(z), 
defined as the n-th derivative of the digamma function 
^o(^) = r'(2;)/r(z). Explicit A-dependencies can be re- 
moved from the right hand side of the fiow equations as 
usual by introducing rescaled dimensionless variables 

. |q| 2nT ~ 5^ ^ T ^ 

One then obtains 

AdJ ^-2S- ^^^^ - || dq g2-3 [,^^1 q, f)] , 

(7) 



AOau = —2u- 



dqif^-^^,[h{5,q,T)], 



(1 + 5f T3 7o 



(8) 

where h{S,q,f) = 1 (1 5)q^-'^/f. The first term in 
each equation is due to the factor A^'^ in the definition of 
the dimensionless variables, the second one captures clas- 
sical fluctuations {ujn = 0), and the third one quantum 
fluctuations (a;„ ^ 0). 



III. SOLUTION OF FLOW EQUATIONS 

For sufficiently small but finite temperature the flow 
passes through two distinct regimes, which are distin- 
guished by the size of the rescaled temperature T. Ini- 
tially one has T <C 1, such that quantum fluctuations 
dominate, while in the final stage, for T ^ 1, the fiow 
is governed by classical fluctuations. In the latter regime 
the third term on the right hand side of the flow equa- 
tions (O and ([H]) can be neglected. For T <^ 1 one can 
use the expansion of the polygamma functions for large 



arguments, ^i{z) 



and ^2{z) ~ — z , to approxi- 



mate the integrals in Eqs. ([7]) and ([8]) as 



(9) 
(10) 



The cutoff scale corresponding to T = 1 is given by Ai = 
(27rT)^/^. Following Millis^ we approximate the flow by 
its quantum contribution with the expansion Eqs. ([9]) and 
pop for A > Ai, and we discard the quantum terms for 
A< Ai. 

The flow equations ([7]) and ([5]) exhibit a fixed point at 
5* = -i, u* = and T* = oo. This fixed point is 
approached by the fiow at the finite temperature phase 
transition, and it describes (approximately) the classi- 
cal non-Gaussian critical fluctuations at Tc- There \5\ is 
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small during the entire flow. Above Tc it remains small 
until the fluctuation contributions to the flow saturate at 
small A. We therefore approximate 1 + (5 ss 1 in the de- 
nominators of the flow equations, which allows us to solve 
them analytically. Note that within this approximation 
fluctuation contributions are symmetric under 5^-^—5, 
while in the exact flow they are larger for (5 < compared 
to (5 > 0, which leads to a suppression of Tc- 

In the quantum regime (A > Ai) the approximate flow 
equations have the form 



-26 , 

zT 

-2u H ^-u"' 

zT 



(11) 
(12) 



Recall that T = 2ttT/A^ is also a flowing quantity. The 
explicit solution for z—3 reads 



vrT 



Cs-A ' 



A'd^ = C^ + /(C3,A), 



(13) 
(14) 



At T=0 one has Ai=0 and the flow can be obtained by 
taking the zero temperature limit of the solution in the 
quantum regime. For z — 3 the unsealed quartic coupling 



saturates at the finite value 

Uq 



Mo = 



1 



1 

2tP 



WqAq 



C3 



(22) 



for A ^> 0. Note that the rescaled variable u vanishes at 
T = 0. For a generic choice of Sq the inverse susceptibility 
scales to a finite value near ^o- At the quantum critical 
point, 



^0 



"0 



/( 



27r2 



,Ao =- 



^^0 



IStt 



2"0 



+ (23) 



the inverse susceptibility scales to zero for A — !> 0. 

For z ~ 2 the quartic coupling vanishes logarith- 
mically for A 0. The inverse susceptibility remains 
generically finite, except at the quantum critical point 
given by 



A2 
^^0 



■ exp 



87r2 
3mo 



Ei 



87r2 
3uo 



(24) 



where 



/(x,A) = iA2-f ixA + ix^ln fl-- 
6 3 6 \ X 



(15) 



The integration constants C3 and C3 are determined by 
the initial conditions at Aq as C3 = ttTAq^Uq^ + Aq and 
C3 = (5o — /(C'sj Aq). The solution for z = 2 is given by 



A^ii^ 



C2 -InA ' 



(16) 



1 



A'^S^ = C2 + ^e2<^^Ei(21nA-2C2) , (17) 



where Ei(x) = J^^dtt^^e* . The integration con- 
stants C2 and C2 are determined by the initial con- 
ditions: C2 = IvrTAo^Uo ^ + InAo and C2 = Sq - 

ie2C=^Ei(21nAo-2C2). 

In the classical regime (A < Ai) the approximate flow 
equations read 



,-,2 



AOaS = -25 - -u 
2 

AOau ~ —2u + 3u 
The explicit solution has the form 
1 



CA2 



A^^^^C'-^ln(A^ + ± 



(18) 
(19) 

(20) 
(21) 



The integration constants C and C" are determined by 
the boundary conditions u^^ = ui and S'^^ — 5i at the 
scale Ai, yielding C = Aj^^ {u'^^ - |) and C = A^^i + 



4C 



2C) 



IV. RESULTS FOR Tc AND Tq 

The phase transition line in the (Jq , T) phase diagram 
is determined by the condition 5^ ^ Q for A 0. Using 
the solution for 5^ in the classical regime, Eq. (|2T|) . this 
yields a condition on the integration constants C and C", 
namely C" — 4^ In ^ . The constants C and C can be 
expressed in terms of the bare variables 5^ and uq by 
matching the initial condition for the classical flow at Ai 
to the solution of the flow in the quantum regime. 

For z = 3, one obtains 



In 



1 



(27rT)i C(T) 



where CiT) = 
temperatures T yields 



27r 
u'T 



2(27rT)2/3 ■ 



4C(T) 

(25) 

Expanding for small 



247r 



T 97r^ 



O 



InT 



(26) 



with Tq = i^f^TT- Note that dependencies on the ultra- 
violet cutoff Ao are absorbed in u'q and Sq'^. Inverting 
Eq. ([25)) to leading order in T yields 

247r ST~Sa 



with ^0 = ^n. 

For 2: = 2, we find 



(27) 



"0 "0 



1 



e^'^^Ei(ln(27rr) - 2C2 



In 



1 + ^dT) 



4C(T) 



(28) 
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where C(T) 



3 



ln(ri/T) with Ti = e^^^/{2Tre 



comes 



(27re) ^AgCxp^^l^y For low temperatures this be- 

0{T/{\nTf) 



T 



7T T\n[\n{T,/T)] 27r 

° ° 3 HTi/T) 3 hi(Ti/T) 



(29) 



Note that the dependence on uq enters only via the scale 
Ti. Inverting Eq. ^ for small T yields 



(30) 



with Ax = f Ti. 

The functional form of the leading temperature de- 
pendencies in Eqs. (|26)) and p9|) is consistent with a nu- 
merical solution of the full flow equations ([7]) and ([S]). 
However, the prefactors are not reproduced exactly, since 
they are affected by the approximations required for the 
sake of an analytic solution. In Fig. 1 we show results 
for the critical line as obtained from Eq. (|25|) for z = 3 
and Eq. (pS)) for z = 2 for specific choice of parameters. 
For z = 3 the logarithmic correction can be clearly seen, 
whereas for z — 2 the (partially compensating) logarith- 
mic corrections are hardly visible, such that the critical 
line looks almost linear. 
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FIG. 1: (Color online): Critical temperature Tc and Ginzburg 
temperature Tg as a function of <5o — ^o"^ for z = 3 (a) and 
z = 2 (b). The bare </)*-coupling is uq — 2, and the ultraviolet 
cutoff Ao = 1 in all cases. The parameter specifying the 
Ginzburg criterion has been chosen as a; = 2/3. The Ginzburg 
temperature below Tc has not been computed here, but has 
been added schematically for completeness. 

We now compute the Ginzburg line in the phase di- 
agram, which marks the boundary of the non-Gaussian 
classical fluctuation regime above Tc- Scaling in the crit- 
ical regime is governed by the interacting Wilson-Fisher 
fixed point (b* , u*). We determine the Ginzburg line from 
the condition {((Aq) = xu* with < a; < 1, where Ag 



is the scale at which 5 crosses zero, that is, 5(Ag) = 0. 
The parameter a: is a measure for the closeness of u to 
the fixed point. This criterion comes about as follows. 
Close to criticality 8 is negative at the beginning of the 
flow, but for r > Tc it eventually increases and diverges 
with the trivial scaling dimension (see Eq. ([rs]) ). If (5 
crosses zero whilst u is still small, u has little influence 
on the flow of 8 also for A < Aq such that we are in 
the Gaussian regime with mean-field exponents. On the 
other hand, if u has reached a value close to its fixed point 
at A = Aq , it affects the flow of b substantially leading 
to non-Gaussian scaling. There is no unique choice of x 
quantifying the "closeness" to the fixed point. This re- 
flects the fact that the Ginzburg line marks a crossover 
regime and not a sharp transition. 

The leading low-T behavior of — computed from 
the Ginzburg criterion described above turns out to be 
the same as that for (5g — S^^ ^ with the same prefactor, ir- 
respective of the choice of x. However, differences appear 
in the first subleading term. For z = 3, one obtains 



and for z = 2, 



-gln(.)T, 



(31) 



(32) 



at low temperatures. We have used the fixed point value 
u* =2/3 as deduced from Eq. (fT9|). Note that the terms 
on the right hand sides are positive. Solving for Tq — Tc 
as a function of b'^ — Jq, one finds that (Tq — Tc)/Tc is of 
order | ln(5g'^ — 5o)|~^ for z = 3, and of order [In | ln((5Q'^ — 
(5o)|]^^ for 2 = 2. Hence, the size of the Ginzburg region 
Tq — Tc is practically of the order near the quantum 
critical point. By contrast, in three dimensions it is of 
order T^i^ The results for Tq are plotted in Fig. [T] for 
the same choice of parameters as above. We can see a 
substantial Ginzburg regime opening between Tq and T^. 
The Tc- and Tc-lines merge when the quantum critical 
point is approached, 5^ 5^^^ since that critical point is 
Gaussian. 



V. CONCLUSIONS 

We have derived analytic expressions for the transi- 
tion temperature Tc and the Ginzburg temperature Tq 
above Tc as a function of the non-thermal control pa- 
rameter (5o near a quantum critical point with a scalar 
(Ising universality class) order parameter in a two dimen- 
sional metal. The calculations are based on fiow equa- 
tions derived from a perturbative renormalization group 
for the Hertz model. The renormalization of the quar- 
tic coupling is crucial to avoid an artificial suppression 
of Tc to zero in two dimensions. Both Tc and Tq are 
essentially proportional to 5q — Jg'^, with logarithmic cor- 
rections depending on the dynamical exponent z. For 
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Tg we confirm the results by Millis.— For Tc we obtain 
the same logarithmic corrections as for Tq , in agreement 
with earlier evidence from a numerical solution of flow 
equations for the symmetry broken phasei^ Neverthe- 
less, the size of the Ginzburg region Tq — Tc, which has 
been calculated analytically for the first time in this pa- 
per, is practically proportional to the distance to the the 
quantum critical point, 60 — 6q'^. Hence, the Ginzburg 
region with its large non-Gaussian classical fluctuations 
covers a substantial part of the phase diagram near a 
continuous quantum phase transition in two dimensional 
metals. Electronic excitations are strongly scattered by 
order parameter fluctuations in that region, which can 



lead to enhanced decay rates, pseudogaps, and other un- 
conventional electronic properties. 
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